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We investigate sequential tunneling through a multilevel quantum dot confining multiple electrons, 
in the regime where several channels are available for transport within the bias window. By analyzing 
solutions to the master equations of the reduced density matrix, we give general conditions on when 
the presence of a second transport channel in the bias window quenches transport through the 
quantum dot. These conditions are in terms of distinct tunneling anisotropies which may aid in 
explaining the occurrence of negative differential conductance in quantum dots in the nonlinear 
regime. 

PACS numbers: 73.21.La, 73.23.-b, 73.23.Hk 



The properties of electron transport in nanoscopic semiconductor devices are an important diagnostic tool in the 
basic study of charge and spin degrees of freedom of functional and reliable nanoelectronic devices [1^ . In quantum 
dots defined electrostatically in a two dimensional electron gas [2], both the electron number and the tunnel coupling 
to the leads are tunable. When high potential barriers separate the quantum dot from the source and drain reservoirs, 
tunneling to and from the dot is weak and the number of electrons N on the dot is a well-defined integer. In 
the Coulomb blockade (CB) regime [3], transport proceeds primarily by sequential resonant tunneling events and 
N fiuctuates by one. In the non-equilibrium regime, where more than one transport channel lies within the bias 
window, a classical transport treatment 4 is no longer appropriate, and a nonequilibrium quantum approach is 
imperative [Sj H]. In this paper, we microscopically derive analytic expressions for the current in the sequential 
regime through the formalism of the reduced density matrix and generalized master equations for the evolution of the 
occupation probabilities [7] . One of the main motivations is the recent observation of negative differential conductance 
in theory and experiment [SJ IHJ [101 IHl HI] in several specific systems. Our goal in the present work is to derive rather 
general but well-defined conditions for the suppression or even the quenching of current whenever a second channel 
enters the bias window. 



We define our model as a two-dimensional lateral quantum dot |13[ll4j weakly coupled to source and drain reservoirs, 
where the total Hamiltonian is given by 



The source and drain Hamiltonians, Hs and respectively, are taken to be noninteracting Fermion systems shifted 
by the bias: 



I. INTRODUCTION 



II. MODEL 



H — Hs + Hqei + Hd + Ht- 



(la) 




(lb) 



with dg,j^ and d^i^j^ a creation and an anihilation operator, respectively, for particles in the source (drain) reservoir. 



The quantum dot Hamiltonian in Eq.(lal is given by. 




(Ic) 
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where the single-particle energies huji are shifted by the applied gate voltage Vg, and Vint is the interaction among the 
confined particles, which we leave general in the present work. 

Finally, the coupling between the quantum dot and reservoirs is described by the tunneling Hamiltonian, 



E 



E 



cI + h.c. 



(Id) 



where h.c. denotes the hermitian conjugate, and where T^^'^^ is the tunneling coeflicient for a particle tunneling 
through the barrier between the single-particle states \s) {\d)) in the source (drain) reservoir and |fc) in the dot. We 
assume operators in each subsystem are independent of each other, and that the eigenstates of the system are known. 
For V\nt 7^ in (Icl, these states are, in general, correlated states — coherent superpositions of Slater determinants. 



III. QUANTUM MASTER EQUATION 



We derive master equations for a non-equilibrium system {Vsd > 0) by means of the reduced density matrix, p„ 
of the system [7j . In the Born-Redfield theory [15j il6j , the diagonal elements are given by 



WmnPnn " E 



nmHmm • 



(2) 



The general form of the transition rate Wmn from state \n) to |m), where both states axe in general many-body 
correlated states, is given by [Tf] . 



= Y E 

a/3 



duj 



~All(Lo, ep - v-)Blp{io + V-) + A±{u, - v+)B^p{u: + v+) 



(3) 



where — e{Vg ± Vb/2), and where 



A±{io,n) = Alt{u:)Q{uj + A^„'l(a;)e(J7 - ^), 
KLi.^) = {n\c^a\'rn){m\cf3\n)5{uj - u)nrn) 



(4a) 
(4b) 



is a generalized quantum dot spectral function with LOnm — ^n^^m and Q{uj) a step function. The analogous spectral 
function for the source reservoir at zero temperature is 



sf/3M = EW'^(^-e^)' 



(5) 



and (w) 

The first term in the integrand of ([3| describes the dot interaction with the source reservoir. This portion of Wn 
may be written as -I- respectively denoting the addition and removal of an electron to or from the dot, 

corresponding to the two terms of ( |4a[ ). The second term in ([S]), describes analogous interactions with the drain 
reservoir, allowing us to symbolically decompose (pi into four terms. 
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denoting the addition and removal of particles from the source and drain reservoirs. 

It is useful to define a generalized chemical potential as the energy required to add a particle to the A^-particle 
ground state, yielding an {N + l)-particle system in the n}^ excited state (n = is therefore the usual chemical 
potential). That is, /ij^ = E'^_^t^ — E^, where is the n}^ excited state of the {N + l)-particle system, and is 

the ground state of the A^-particle system. At zero temperature and zero bias, the system is in the iV-particle ground 
state when the chemical potential /ijv-i '^^'^^ below the Fermi energy of the reservoirs, and /ijj^ lies above the Fermi 
energy. 



IV. TRANSITION CONDITIONS 



If we denote by N the number of confined particles at equilibrium, and by N^. the number of particles in the state 
|fc), then by an analysis of Eqns. (4a I, (4b) and ([5|, one can deduce that Wmn is non-zero only when one or more of 
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the following four sets of conditions are satisfied: 



N<N„, 




+ 1, (4"^ 


-3+ 


) < (A„,„± ieFs), 
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-) > (A„„ T ^e^s), 


(7b) 


N N„, 
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-) > (A„„± ieFs), 


(7c) 


N < iV,„ 




-1' K:'- 


^9 


< (A„„ T ^eT^s). 


(7d) 



Here, i^ff = f/^j — ^% denotes the difi^erence in chemical potentials between M particles and N particles in the dot. 
In addition, = — {ep — eVg) denotes the energy difi'erence between the chemical potential of the iV-particle 
dot and the chemical potential of the source reservoir, and similarly g~ = (ep — sVg) — fJ,%_i denotes the energy 
difference between the chemical potential of the drain reservoir and that of the {N — l)-particle dot. The energy 
difference Amn = SE"^ — 5E"-, where (5i?" is the excitation energy of the state \n). That is, if is the ground- 
state energy of the A^fc -particle system, then E^ = Eg" + SE*" with SE'' ^ 0. The relations iV ^ iV^ = iV„ ± 1 
in Eq. ([7) are a consequence of sequential tunneling by (only) a single electron. On the other hand, the relations 
i'^Ni + y ) i^pq i eVB/2) are a consequence of resonant tunneling. Moving out of the Coulomb blockade regime 
relaxes the first set of inequalities while non-Mar kovian effects relax the second set. 



V. CURRENTS: 1 AND 2 CHANNEL SYSTEMS 



We consider the current through a quantum dot such that electrons tunnel into the dot from the source reservoir, 
and out of the dot to the drain. The evolution of the diagonal elements of the density matrix (|2| for the system 
depends on both the source and drain reservoirs and can be separated into a contribution due to coupling with the 
source, and a contribution due to coupling with the drain, 

Pmm = Pmm ~^ Pmm ~ ^ ^ [(^mraP"" ~ ^nmPnim) + {Wmn Pnn ~ ^nm Pmm)] • (8) 
n 

We denote current into the dot through the source barrier by 7'^, and current out of the dot through the drain 
barrier by I^. At long times, the current is proportional to the difference between these two quantities: I = — = 
e{N)s - e{N)D, where, for example, {N)s = E,„^mPm- 

Using ([s]) we obtain, 

= e ^ ^rn {W^tPn " W^+Pm) (9a) 

mn 

I^ = eY,Nm {W^^Pra - W^;;pn) , (9b) 



mn 



and = — in the steady-state regime. 

We consider the configuration shown in Fig. [T] where the initial equilibrium state of the system is the iV-particle 
ground state, and we calculate the steady state currents through the dot for the respective cases of one (Fig. [T|d) 
and two (Fig.flj:) channels in the transport window. We define three many-body states: |1) = \N)o is the A^-particle 
ground state, \2) = \N + l)o is the {N + l)-particle ground state, and |3) = |7V)i is the A^-particle first excited state. 
The transitions between these states are governed in part by Eq. ( [7| . 

There will be precisely one channel in the transport window, Fig.TTlD, whenever all of the following three inequalities 
are satisfied: 

g+ < ^eVB < g- (10a) 
^eVB<A^-g+ (10b) 
^eFs < A^+i+5+, (10c) 
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where A^ = E^^t^,^,. — -E'^ound excitation energy of the first excited state of the iV-particle system (this 

excitation can be due to an electron-electron interaction, or due to Zemann spliting in the presence of a magnetic 
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FIG. 1: Diagram of dot energy regimes, (a) No transport channels in the bias window, (b) single transport channel, p5v, in the 
bias window, (c) two transport channels, /i'^ and /xjv-i, in the bias window. 



field [18j). Assuming initially unit occupation of only the A^-particle ground state, the infinite coupled set (|2| reduces, 
on application of Q to two equations: 



Pii 
P22 



W12P22 
W21P11 



W21P11, 

Wi2P22- 



(Ua) 
(lib) 



That is, transport involves the ground states of the N and {N + l)-particle systems (only), as expected. Solving this 
coupled set and inserting into ([9]), the steady state currents are found to be 



s+ 



(12) 



which is the standard expression 1191 EU] for a two- level steady state dc current. 

The two-channel regime. Fig. [ip, is entered upon relaxing condition (10b I. If we again assume unit occupation of 
the A^-particle ground state as our initial condition, the set Q now reduces to three equations: 

Pii = 1^12^22 - W21P11 (13a) 

P22 = W21P11 + VP23P33 - P22{Wi2 + ^^32) (13b) 
/O33 = M^32P22 - VF23P33 (13c) 

Transport now involves the ground states of the N and (Af-l-l)-particle ground states — states |1) and |2) respectively — 
as well as the first excited state with TV particles — state |3). (Note that if condition ( |10c[ ) rather that (10b I were 
relaxed, the expressions, and the states involved, would differ.) Solving this coupled set and inserting into (|9|, the 
steady state currents are found to be 



= If = I? = 



S+i 



(14) 



Upon increasing the bias from zero, the steady-state transport current should change from zero to lb (12), and 



finally to /c (14 1. Negative differential conductance will be observed whenever Ic < h- That is, an increase in the bias 



reduces the current. General conditions under which such phenomena occur are explored in the subsequent sections. 



VI. ASYMMETRIES 



To explore the relative magnitudes of and 1^ we define two anisotropies. The first is an extrinsic anisotropy 
characterizing the barrier widths. If we assume tunneling to and from the source is proportional to tunneling to and 
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from the drain, we can write 

W^tlW:^,n ^ «mn (15) 

for all many-body states \n) and \ra). Through Eq. ([s]), these can be related directly to the tunneling coefficients T^, 
Ti appearing in Eq. (Id). In particular, the noninteracting limit yields 



limit 



where and T'^° are the tunneling amplitudes of the resonant states. In general, < a < oo. 

The second anisotropy is intrinsic to the dot and denotes the relative ease with which particles can tunnel to 
different orbitals. Core orbitals tightly bound to the center of the dot, for example, have weaker coupling to the leads 
than do orbitals localized around the edge of the dot. for the particular case of the many-body states |1), |2), and |3) 
defined above, we write 

W^'^-IW^^' EE e, (17) 

where |1) and |3) are A^-body states, and |2) is an {N + l)-body state. The magnitude of e will be determined by the 
spatial extent of the single-particle orbitals involved in the respective transitions. Starting from the {N + l)-particle 
state |2), describes an electron tunneling out of the dot, leaving the system in the A^-particle state |1), whereas 

W^2~ describes a similar process, with the system ending up in state |3). The ration e depends on the particular 
orbitals involved. 

In the example we give below, 1 < e < oo. 



VII. QUENCHING OF CURRENT 



Examining the currents Ic, Eq. (|14|, and 7^, Eq. (|12|, in terms of the anisotropy parameters a- 

= R. 



Ic 

h 



(l + e)(l-|-a2i) 



e 1 -I- Q!2i 



021 
"23 



mn and e, we find 
(18) 



For i? < 1, current is suppressed and, in the extreme case of q;23 — > (i.e., IV^^ 0), the presence of the second 
channel in the transport window completely quenches the current, /c ^ 0. A detail of the Coulomb diamond for 
which this analysis applies is given in Fig. [2] 

As an illustrative example, we consider a two-dimensional parabolic quantum dot, confining two electrons, with se- 
quential resonant tunneling proceeding by the addition and removal of a third electron. The three states corresponding 




FIG. 2: Coulomb diamond detail showing suppression of ground-state current, the current increases along (i), but may decrease 
along iii). In extreme cases, the decrease along {ii) matches the increase along (i), yielding no current along [Hi). It and Ic 
correspond respectively to (b) and (c) of Fig. [l] and their analytic expressions are given in Eqns. ( |12[ ) and ( |14[ | respectively. 




FIG. 3: Configurations capable of quencliing ground-state transport. Rates involving states |1) and |2) involve tunneling of an 
electron into a dot edge-state (more probable) while those involving states |3) and |2) involve tunneling into a dot core-state 
(less probable). 



to Fig.[T|: are the ground-state singlet, |1), the lowest-energy triplet, |3), and the three-particle ground state, |2). These 
respective configurations are shown in Fig.|3] The orbitals can be labeled by two harmonic oscillator indices plus spin, 
\mns), where the shell number is given bym-|-n = 0,l,2,..., and the angular momentum = n — m = 0,±1, ±2, .... 
With regard to Fig. |3) we have |1) = |00 T,00 i), |2) = |00 T,00 i,01 i), and |3) = |00 i,01 i). The orbital |01) 
penetrates the barrier [T] to a greater degree than the |00) orbital, and so the tunneling amplitude Tqq'^^ has smaller 
magnitude than Tqi"^"^; depending on the geometry and the magnetic field, this difference can be severe. In general. 



we expect e > 1 in Eq. (17 1 



In the two-channel regime, we look first at the symmetric barrier case, amn = 1- Here, suppression of current, 
R < 1, Eq. (18 1, occurs whenever e > 2, that is, whenever tunneling into the edge is more than two times as likely as 
tunneling into the core. In this case, current is never completely quenched; the limiting case (e —^ oo) has R — > 2/3. 

Complete quenching does occur, however, if we additionally have a barrier asymmetry. Specifically, if the source 
barrier is sufficiently opaque such that tunneling into the core from the source reservoir is suppressed (W;^"*" — > 0, see 
Fig. [3 1, but tunneling from the core to the drain is not {W^~ ^ 0), then 023 and /c ^ even though If, ^ 0; 
the addition of the second transport channel in the bias window suppresses the current completely. Once state |3) 
becomes occupied, the system is blocked and no other transitions can occur. 

This quenching of current is a dynamical blockade: it requires an interplay between both ground and excited states. 
This is in contrast to effects such as a) spin blockade, where a transition changing the number of electrons by one is 
accompanied by a change of spin greater than 1/2 |21j . or b) blockade due to spin polarized injection and extraction 
of charge carriers |22j . where the current is quenched due to the Pauli exclusion principle. 

The quenching of current will be destroyed whenever the spin-flip plus orbital decay rate is large compared to 
the dwell time of the tunneling electron. In this case, a direct transition from |3) to |1) may occur by emission of 
a phonon and simultaneous absorption of spin through nuclear baths, spin-orbit coupling, or other spin-symmetry- 
breaking process. Likewise, strong Coulomb correlations will generally open up additional tunneling pathways, likely 
involving edge orbitals, and this too will destroy the quenching of current. 



VIII. CONCLUSIONS 



In conclusion, we have derived a set of general yet well-defined conditions for quenching of transport through a 
multichannel quantum dot. We find that under certain tunneling anisotropics the presence of a second transport 
channel may decrease the current through the dot. Specifically, if the tunneling rates to or from core states are less 
than half than the rates to or from edge states, the current may be suppressed. Complete quenching of the current 
may occur if, in addition to the edge/core anisotropy, an anisotropy between source and drain barriers is also present. 
The presence of these dark channels is a measure of the lifetime of excited states. 

Finally, although the present work focused on the two-electron case, the general phenomena should appear in 
numerous locations throughout the Coulomb diamonds; it can generally appear whenever ground-state tunneling 
transitions involve edge orbitals, and excited-states tunneling involves core orbitals, and asymmetric source and drain 
barriers are present. 
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